15/WO 73 


NASA TECHNICAL NOTE 


tn 

oo 

r>^. 



<1 

CO 


■< 


DESIGN OF EXPERIMENTS FOR MEASURING 
HEAT-TRANSFER COEFFICIENTS WITH 
A LUMPED-PARAMETER CALORIMETER 

G. James Van Fossen, Jr. 

Lewis Research Center and 
U.S. Army Air Mobility R&D Laboratory 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D 




f 

f 

NASA TN D-7857 





C. • JANUARY 1975 



1. Report No. 2. Government Accession No. 

MSA TN D-7857 


4. Title and Subtitle 

DESIGN OF EXPERIMENTS FOR MEASURING HEAT -TRANSFER 
COEFFICIENTS WITH A LUMPED -PARAMETER 
CALORIMETER 


7. Author(s) 

G. James Van Fossen, Jr. 


9. Performing Organization Name and Address 

NASA Lewis Research Center and 

U.S. Army Air Mobility R & D Laboratory 

Cleveland, Ohio 44135 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


3. Recipient's Catalog No. 


5. Report Date 

January 1975 


6. Performing Organization Code 


8. Performing Organization Report No. 

E-8104 


10. Work Unit No. 

505-04 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Note 


14. Sponsoring Agency Code 



16. Abstract 

A theoretical investigation was conducted to determine optimum experimental conditions for 
using a lumped -parameter calorimeter to measure heat-transfer coefficients and heating rates. 
A mathematical model of the transient temperature response of the calorimeter was used with 

•J 

the measured temperature response to predict the heat-transfer coefficient and the rate of 
heating. A sensitivity analysis was used to determine the optimum transient experiment for 
simultaneously measuring the heat addition during heating and the convective heat-transfer 
coefficient during heating and cooling of a lumped-parameter calorimeter. Optimum experi- 
ments were also designed for measuring the convective heat-transfer coefficient during both 
heating and cooling and cooling only. 


17. Key Words (Suggested by Author(s) ) 

Heat-transfer coefficient 
Sensitivity analysis 
Optimization 


18. Distribution Statement 

Unclassified - unlimited 
STAR category 34 (rev. ) 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

19 


22. Price* 
$3.00 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 
























DESIGN OF EXPERIMENTS FOR MEASURING HEAT -TRANSFER COEFFICIENTS 
WITH A LUMPED-PARAMETER CALORIMETER 
byG. James Van Fossen, Jr. 

Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

A theoretical investigation was conducted to determine optimum experimental cbn- 
ditions for using a lumped-parameter calorimeter to measure heat-transfer coefficients 
and heating rates. A mathematical model of the transient temperature response of the 
calorimeter was used with the measured temperature response to predict the heat- 
transfer coefficient and the rate of heating. 

A sensitivity analysis was used to find optimum experimental conditions for the 
lumped -parameter calorimeter. An optimum transient experiment was designed for 
simultaneously measuring the heat addition during heating and the convective heat- 
transfer coefficient diming heating and cooling of the calorimeter for steady convective 
-conditions. 

Optimum experiments were also designed for measuring the convective heat- 
transfer coefficient during both heating and cooling and cooling only. 

For the case of measuring the heat -transfer coefficient from the cooling of the 
calorimeter only, the optimum experiment consists of allowing the dimensionless tem- 
perature of the calorimeter to fall from 1.0 to about 0.183. For this case the rate of 
heat addition to the calorimeter during heating has no effect on the optimum experiment. 

For the other two cases, that is, the cases where both heating and cooling data of 
the calorimeter were used, it was shown that the optimum experiment was dependent on 
dimensionless heating rate. A heating rate as slow as practical was shown to be the 
best. The optimum dimensionless duration of the experiment was shown to be dependent 
on heating rate for the case of measuring both heat-transfer coefficient and heating rate 
simultaneously and the case of measuring heat-transfer coefficient only . 



INTRODUCTION 


An analytical investigation based on sensitivity coefficients was conducted to deter- 
mine the optimum transient experiment for measuring convective heat-transfer coef- 
ficients with a lumped -parameter calorimeter (LPC). Lumped-parameter calori- 
meters consist of a high -thermal -conductivity material having a known specific heat. 
They can be imbedded in the surface of a larger body to measure local heat-transfer 
coefficients. The usual approach for these transient tests is to have the gas stream 
conditions at steady state and to heat the calorimeter electrically to a predetermined 
temperature above that of the surrounding surface. The heater is then switched off, 
and the calorimeter is allowed to cool by convection. By measuring the temperature of 
the calorimeter as a function of time during the heating and cooling cycles, it is possible 
to determine, by use of a mathematical model, the heat added by the electric heater 
and the heat-transfer coefficient on the surface at the location of the LPC. 

Heat-transfer investigations that have used the LPC approach are reported in ref- 
erences 1 to 3. In all three references, only the cooling portions of the cycle were used 
for determining convective heat-transfer coefficients. In addition, the coefficients in 
references 1 and 2 were based upon only two temperature measurements - the initial 
temperature of the LPC prior to beginning the cooling portion of the cycle, and a single 
time -temperature measurement during the cooling cycle. The coefficients obtained in 
reference 3 were based upon the slope of the experimentally measured time -temperature 
curve. 

The present investigation was conducted to show the optimum time for making tem- 
perature measurements and to introduce the method of least squares for utilizing multi- 
ple time -temperature measurements during both the heating and cooling cycles for in- 
creased accuracy in heat -transfer -coefficient and heat -addition measurements. 

Optimum experiments are defined for the following three cases: (1) heat -transfer co- 
efficients calculated by using cooling data only, (2) heat -transfer coefficients calculated 
by using both heating and cooling data, and (3) both heat -transfer coefficients and heat 
addition to the calorimeter during the heating cycle calculated simultaneously from 
heating and cooling data. 


THEORY 

In this section a mathematical model is developed for the LPC. It is also shown 
that the least-square method can be used to calculate the convective heat-transfer coef- 
ficient from the data and the mathematical model. 

The LPC is usually made of a high -thermal -conductivity material such as copper 
or aluminum. The high thermal conductivity helps to promote a uniform temperature 
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distribution. The assumption of uniform temperature allows the calorimeter to be 
lumped into a single mass. This simplifies the mathematical modeling of the calori- 
meter . 

In order to maintain a uniform temperature in the calorimeter, the Biot number 
must be kept less than about 0.1 (ref. 4). The Biot number is given by 

Bi = ^ (1) 

k 

(All symbols are defined in the appendix.) 

Equation (1) shows the importance of using a high -thermal -conductivity material for 
the calorimeter. The characteristic length L is obtained by dividing the volume of the 
calorimeter by the convection surface area. Thus, the Biot number can also be c<?n- 
trolled by the physical dimensions of the calorimeter. 

The LPC has some means to input heat, usually an electric heater. A thermocouple 
is usually used to measure the temperature of the calorimeter. The operation of the 
LPC is quite simple. The convective flow is established, the heater is turned on, and 
the calorimeter is heated to some temperature above ambient. When the predetermined 
maximum temperature is reached, the heater is switched off and the LPC is allowed to 
cool by convection. The temperature is recorded as a function of time during both heat- 
ing and cooling. A mathematical model is then used to estimate the convective heat- 
transfer coefficient. 


Mathematical Model j 

i *j 

The assumptions used for a mathematical model of the LPC are 

(1) Temperature in the calorimeter material is uniform. 

(2) Heat -transfer coefficient remains constant throughout both heating and cooling 
of the LPC . 

(3) Calorimeter material properties are known. 

(4) Ambient temperature remains constant. 

(5) Air properties remain constant. 

(6) Heat addition to calorimeter is constant until heater is switched off. 

With these assumptions a simple heat balance gives the dimensionless temperature his- 
tory of the LPC as 
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where t q is the dimensionless time at which the heater is switched off. From equa- 
tions (2) and (3) t q must be 



The dimensionless temperature $ is defined as 
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where ^T max is the maximum temperature rise above ambient. This temperature ris< 

is determined before the experiment starts and is an experimental constraint; the heater 

is not switched off until AT„„ is attained. 

max 

The dimensionless heating rate /3 is defined as 


Q 


hA AT 


max 


( 6 ) 


Note that if 8 is less than or equal to 1.0, AT will never be attained. The dimen- 

IllaX 

sionless time r is given by 


r CO 

P LC p 

Figure 1 shows dimensionless temperature as a function of dimensionless time for sev- 
eral dimensionless heating rates. 


Calculation of Heat-Transfer Coefficient 

One method that is used to calculate the heat-transfer coefficient from the data is to 
use a single time -temperature measurement during the cooling of the LPC. The 
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temperature of the LPC above ambient is assumed to be known at time r . The experi- 
mentally measured temperature and time at some time greater than t q are inserted 
into equation (3), and the heat -transfer coefficient is determined from the resulting 
equation . 

A more accurate method of calculating h from the LPC would be to use more than 
one data point. The method of least squares allows the use of any number of data points 
in the calculation. The method of least squares consists of minimizing the least-square 
function 

N 

F < h >-£[ T i- T i,ci< h f < 8) 

i=l 

with respect to the unknown heat-transfer coefficient h. The value of h that minimizes 
F(h) is the best value. 

The method of least squares allows all the data to be used. With this method, data 
from the heating portion of the cycle as well as the cooling data can be used. If both the 
heating and cooling data are used, it is also possible to simultaneously calculate the 
heating rate and the heat -transfer coefficient. 


ANALYSIS OF OPTIMUM EXPERIMENT 

An optimum experiment is defined as one which allows the parameters in a math- 
ematical model to be calculated from data with the greatest accuracy. One might think 
that experimentalists always perform experiments that will give the greatest accuracy; 
however, this is not true in the sense meant here. One can run the best experiment in 
terms of careful specimen preparation, placement of sensors, test procedure, and data 
acquisition and still have a poor experiment for measuring parameters. It is demon- 
strated by example later in this section that errors in the calculated heat-transfer co- 
efficients caused by errors in the data can be amplified if measurements are taken at 
the wrong time. The effect of errors in the data can be minimized if measurements are 
taken at the optimum time. 

The experimental conditions which yield predicted temperatures that are most sen- 
sitive to changes in the parameters are the optimum conditions. The experimental con- 
ditions which can be changed (controllable variables) in LPC tests are maximum tem- 
perature rise, characteristic length of the calorimeter, heating rate, duration of the 
experiment, and calorimeter material. 

Sensitivity coefficients can be used to find the optimum experimental conditions. A 
sensitivity coefficient is defined as the derivative of the predicted value of a measured 
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quantity with respect to the parameter to be estimated. Sensitivity coefficients are made 
dimensionless by multiplying by the parameter and dividing by a function of the measured 
quantity. Sensitivity coefficients show how the quantity predicted from the mathematical 
model is affected by small changes in the parameters. 

For the LPC the parameters we wish to estimate are the convective heat-transfer 
coefficient h and the heating rate Q. For the case of predicting h from the cooling 
of the calorimeter, the sensitivity coefficient is obtained from equation (3) by differenti- 
ating with respect to h. The time r Q is set to zero because the heating cycle does not 
affect the temperature for this case. We obtain in dimensionless form 


S T - 
\ - 


d -l=-re~ T 


AT 9h 
max 


(9) 


Figure 2 shows the dimensionless sensitivity coefficient for this case as a function of 
dimensionless time. 

T 

For both heating and cooling the dimensionless sensitivity coefficient is given by 



-/3 (1 - e -T ) + /3Te~ T T — T 0 



( 10 ) 
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T 

Figure 3(a) shows as a function of r. 

For both heating and cooling the dimensionless sensitivity coefficient 



is given by 


( 12 ) 


(13) 
T 

Figure 3(b) shows Sq as a function of r. 

Extending the work of Box and Lucas (ref. 5), Beck (ref. 6) proposed a criterion 
to find the optimum experimental conditions. For a fixed temperature range and a large 
number of equally spaced measurements in time, Beck's criterion is to maximize the 
determinant of the matrix D. The elements of D are given by 


f/3(l - e -T ) t<t q 




( 14 ) 
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where r m is the duration of the experiment. The maximum temperature rise is a con- 
straint and must be reached sometime between zero and t . 

The optimization procedure then is to fix all the controllable variables except r m 
and calculate 

D = det D (15) 

as a function of r m , where D is the sensitivity criterion. The maximum value of D 

max[D] = D opt (16) 

T m 

indicates the best value of T m for that particular setting of the other controllable vari- 
ables. This procedure is repeated for different settings of the other controllable varia- 
bles until the best combination is found. 

As an example of how the sensitivity criterion D is formed, consider the case of 
predicting the two parameters Q and h. Let 

Pj=h (17) 


and 

P 2 = Q (18) 

Then from equation (14) the elements of the matrix D are 
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The sensitivity criterion D is then 

5 = d ll d 22 - < d 12> 2 < 22 > 

In order to demonstrate quantitatively the advantage of the least-square method over 
the single -point method, a sample case was constructed. The case of measuring the 
heat -transfer coefficient from cooling was chosen. The mathematical model for the di- 
mensionless temperature is 


$ = e~ at (23) 

The exact value of the parameter a was taken as 1.0; thus, dimensionless time cor- 
responded to real time. The effect of errors in the measured value of on the calcu- 
lated value of the parameter ct was investigated by using the single -point method to 
calculate a from simulated data. The simulated data were generated from the model 
(eq. (23)) by randomly adding or subtracting an error of 1 percent of the maximum value 
of $ (0.01) to the calculated <£ for each time. The parameter o> was calculated from 

a = _ ln(* ± error) ( 24 ) 

t 

for each simulated data point. 

The least-square method was also used to calculate a from measurements con- 
taining errors. For this case the same 1 percent of maximum $ error was randomly 
added to or subtracted from the exact value of $ for seven data points that were equally 
spaced from 0, 25 to 1. 75 units of time. 


RESULTS 

For the case of calculating h from cooling only, equation (3) shows that the only 
controllable variable which appears in the model is dimensionless time. For this case, 
only a single curve is necessary for the optimization. Figure 4 shows the sensitivity 
criterion as a function of dimensionless duration of experiment. The optimum dimen- 
sionless duration is about 1.7 after the heater is switched off. From the mathematical 
model (eq. (3)) this corresponds to letting the dimensionless temperature fall from 1.0 
to about 0.183. For this case the sensitivity criterion is not affected by either the heat- 
ing rate or the maximum temperature rise. The experimenter should let common sense 
dictate the maximum temperature rise. Too small a temperature rise will allow noise 
and other measurement errors to degrade the data. 
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For the case of calculating h from both the heating and cooling data the sensitivity 
criterion becomes a function of the dimensionless heating rate as well as the dimension- 
less duration of the experiment. Figure 5 shows the sensitivity criterion as a function 
of dimensionless duration of experiment for several values of dimensionless heating 
rate. From this figure it appears that having the dimensionless heating rate near unity 
is a desirable condition. 

If the dimensionless heating rate (3 were unity, the LPC would take an infinite time 

to reach AT . The optimum dimensionless time t , as a function of f3 is shown 
max opt 

in figure 6. The experimenter should choose 0 as near unity as possible, but he must 
maintain a reasonable real time t. The dimensionless time is given by equation (7). 

The real-time duration of the experiment can be shortened for a given h by making 
L small. This is also consistent with the requirement of uniform temperature (eq. (1)). 

The sensitivity criterion is again a function of both heating rate and dimensionless 
time for the case of calculating both Q and h from heating and cooling. Figure 7 
shows the sensitivity criterion as a function of dimensionless time for several heating 
rates. The best experiment is again one in which /3 is near unity. 

Figure 8 shows the optimum dimensionless time as a function of dimensionless heat- 
ing rate for the case of calculating both h and Q from the heating and cooling data. 

The results are similar to the case of calculating h from heating and cooling data. 

The results of the sample case used to demonstrate the advantage of the least-square 
method over the single-point method are shown in figure 9. The upper and lower curves 
in the figure show the errors in the parameter a calculated by using the single -point 
method (eq. (24)) as a function of the time the measurement was taken. For the upper 
curve the error was subtracted from $, and for the lower curve the error was added. 

The data used to calculate the parameter a with least squares are shown in fig- 
ure 10. The circular symbols are the simulated data points, and the curve is the math- 
ematical model of equation (23). The resulting error in a calculated by least squares 
is shown in figure 9 as a constant over the range of measurement. The error is 
0.226 percent. Figure 9 verifies that greater accuracy can be obtained with the least- 
square method than with the single-point method. Figure 9 also shows that with the 
single-point method the effect of measurement error can be minimized by choosing the 
correct time to take the measurement. The figure indicates that the error in the cal- 
culated value of a is minimum when the time is 1.0. This is the same time at which 
the sensitivity coefficient shown in figure 2 is a maximum. 


CONCLUSIONS 

A sensitivity analysis has been used to design the optimum experiment for a lumped- 
parameter calorimeter (LPC) used to measure heat -transfer coefficients. Three cases 
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have been considered. In the first case the LPC was used to measure the heat-transfer 
coefficient from cooling data only. Both the heating and cooling cycles were used to 
measure the heat -transfer coefficient in the second case. For this case the heat input 
during heating was assumed to be known. Finally, for the third case, both the heating 
and cooling data were used to measure both the heat-transfer coefficient and the heat 
input during heating. 

The optimum conditions for the first case were not affected by the heating rate. The 
optimum dimensionless time for this case was about 1.7. This corresponds to letting 
the LPC cool off until the dimensionless temperature is about 0.183. 

For the second case, both the heating and cooling cycles were used to measure the 
heat-transfer coefficient. The sensitivity criterion is affected by the heating rate for 
this case. The dimensionless heating rate should be chosen near unity. It may not be 
practical to choose the dimensionless heating rate too near unity because the real-time 
duration of the experiment may become too large. The real-time duration of the experi- 
ment can be controlled by varying the physical dimensions of the calorimeter. Once the 
heating rate is chosen, the optimum dimensionless duration of the experiment can be 
obtained. The mathematical model can then be used to translate this time into a tem- 
perature . 

For the third case, both the heating and cooling cycles were used to measure the 
constant heat-transfer coefficient and the constant heating rate during heating. The opt- 
imum dimensionless heating rate was found to be near unity, as in the second case. The 
real-time duration of the experiment can be made reasonable by designing the calori- 
meter such that its characteristic length is small. Once the dimensionless heating rate 
is chosen, the optimum dimensionless time can be found. 

A comparison of the single-point and multiple -point (least square) methods of meas- 
uring the heat-transfer coefficient from the cooling of the LPC shows that the least- 
square method results in less error. The error in the heat -transfer coefficient calcu- 
lated by using the single -point method can be minimized if the measurement is taken at 
the time when the sensitivity coefficient is maximum, that is, when dimensionless time 
is 1.0. 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U.S. Army Air Mobility R&D Laboratory, 

Cleveland, Ohio, November 5, 1974, 

505-04. 
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APPENDIX - SYMBOLS 


A 

Bi 

c f 

D 

D 

D 


opt 


d « 

h 

k 

L 

N 

Pi 

Q 


T i 


T ; 


i, cal 


s 

AT 

t 


max 


a 

0 

P 

T 


T. 

1 


T 


m 


T 

O 


area for convection 

Biot number, hL/k 

specific heat of calorimeter material 

sensitivity matrix whose elements are 

determinant of matrix D 

maximum value of D as defined by eq. (16) 

elements of D defined by eq. (12) 

convective heat-transfer coefficient 

thermal conductivity of calorimeter material 

characteristic length of calorimeter material, volume/area 

number of observations 
th 

l parameter in a mathematical model 
heat added to calorimeter per unit time 

sensitivity coefficient, derivative of predicted temperature with respect to 
parameter p^ 

temperature of lumped -parameter calorimeter 
measured temperature at time t. 

temperature calculated from mathematical model at time t. 
ambient temperature 

maximum temperature rise above ambient 
real time 

example time constant, in eq. (23) 
dimensionless heating rate, Q/hA AT 

IT) cLX 

density of calorimeter material 
dimensionless time, ht/pLC 

discrete dimensionless time at which T^ is measured 
dimensionless duration of experiment 
dimensionless time at which heater was switched off 



T opt optimum duration of experiment 

$ dimensionless temperature of lumped -parameter calorimeter (T - T )/AT 

s max 
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(t>) Sensitivity coefficient, Sq, 


Figure \ - Effect of heating rate on sensitivity coefficients for both heating and cooling. 
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Figure 6. - Optimum dimensionless time as (unction of heating rate (or estimating 
heat-transfer coefficients from heating and cooling data. 


Dimensionless 



Figure 7. - Sensitivity criterion as function of time for estimating 
heat-transfer coefficients and heat addition from heating and 
cooling data at several heating rates. 



Figure 8. - Optimum time as function of heating rate for estimating both heat -transfer 
coefficients and heat addition from heating and cooling data. 
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Dimensionless temperature, $ Error in example time constant a, percent 



Figure 9. - Error in estimating example time 
constant with a 1 percent measurement 
error as function of time at which meas- 
urement is taken. 



Figure 10. - Simulated data used in least- 
square example. 


NASA- Langley, 1975 


E-8104 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

I BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 




POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 





